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Abstract 

We consider the ground state energy of a spinor field in the background 
of a square well shaped magnetic flux tube. We use the zeta- function regu- 
larization and express the ground state energy as an integral involving the 
Jost function of a two dimensional scattering problem. We perform the 
renormalization by subtracting the contributions from first several heat 
kernel coefficients. The ground state energy is presented as a convergent 
expression suited for numerical evaluation. We discuss corresponding nu- 
merical calculations. Using the uniform asymptotic expansion of the spe- 
cial functions entering the Jost function we are able to calculate higher 
order heat kernel coefficients. 



Introduction 



Since the classical work of H.B.G. Casimir, where the energy of vacuum polar- 
ized by two conducting planes was calculated, a number of similar problems was 
investigated for various external conditions (boundary conditions, background 
potentials etc.), see, for example, the recent review p| or the books 0, Q . No 
general rule for the dependence of the vacuum energy on the background proper- 
ties has been found so far. In particular, it is unknown how to forecast the sign 
of the energy. 

An interesting problem is the calculation of the vacuum energy of a spinor 
field in the background of a magnetic field. This problem was first considered as 
early as in 1936 by W. Heisenberg and H. Euler ^ who were interested in the 
effective action in the background of a homogeneous magnetic field. 
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Presently, the interest is shifted to string hke configurations. In the con- 
tribution of the fermionic ground state energy to the stabihty of electro weak 
strings was addressed. In |^ the gluonic ground state energy in the background 
of a center-of-group vortex in QCD had been considered. 

A very special inhomogeneous refiectionless magnetic background is for ex- 
ample the model for domain wall p]. In this case, the equation for the vacuum 
fluctuations was easy to solve analytically. However, the final expression for the 
effective action involves an additional integration over a momentum, thus bring- 
ing no real simplifications. 

In QED, for the background of a fiux tube with constant magnetic field inside 



the ground state energy was calculated in []T3l. Here it turned out to be negative 



remaining of course much smaller than the classical energy of the background 
field. 



An interesting approach is that used in |]TD|] where the issue of the sign and of 
bounds on fermionic determinants in a magnetic background had been considered. 

The remarkably simple case of a magnetic background field concentrated on a 
cylindrical shell (i.e., a delta function shaped profile) was considered in [|ll|. Here 
different signs of the vacuum energy turned out to be possible in dependence on 
the parameters (radius and flux). 

In general, the issue of stability of the strings is of interest. In electro weak 
theory, and in QED in particular, the coupling is small. Hence the vacuum 
energy as being a one loop correction to the classical energy is suppressed by this 
coupling. While in QED a magnetic string is intrinsically unstable, in electro weak 



theory there are unstable and stable conflguration known [12|. Here quantum 



corrections may become important for the stability. A question of special interest 
is whether a strong or singular background may have a quantum vacuum energy 
comparable to the classical one. 

For the calculation of ground state energy, it is necessary to subtract the 
ultraviolet divergences. This procedure is in general well known and we follow 
[]. Roughly speaking one has to subtract the contribution of the flrst few heat 
kernel coefficients {a^ through 02 in the given case of (3+1) dimensions). After 
that one can remove the intermediate regularization and one is left with flnite 
expressions. However, in order to obtain these flnite expressions in a form suitable 
for numerical calculations one has to go one more step. As described in [], see 
also below in section 2, one has to add and to subtract of a certain part of the 
asymptotic expansion of the integrand. 

In the present paper we generalize the analysis done in to a rectangular 
shaped magnetic background fleld. So we consider the vacuum energy of a spinor 



fleld in QED for a background as given by Eqs. (2.1, 2.5). This problem is 



technically more involved and allows progress in two directions. First, it allows 
to reflne the mathematical and numerical tools for such problems and, second, it 
allows to address the question how the vacuum energy behaves for an increasingly 
singular background (making the rectangle narrower). So this model interpolates 
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to some extend between the flux tube with homogeneous field inside in H and 



the delta shaped one in 



The basic principles of this procedure are described briefly in Sec. 2. Namely 
it starts with the well known zeta functional regularization. The regularized 
ground state energy is represented as a zeta function of a hamiltonian spectrum 
and treated in termini of heat kernel expansion. The representation of regularized 
ground state energy as an integral of the logarithmic derivative of the Jost func- 
tion for spinor wave scattering problem on the external magnetic background 
is obtained. The Jost function is obtained from the exact solutions of Dirac 
equation, derived in Sec. 3. The explicit form of the exact and asymptotic Jost 
function is considered in Sec. 4. In Sec. 5 the representation useful for further 
numerical evaluations is derived. Sec. 6 is devoted to the calculation of the heat 
kernel coefficient 05/2 and the Sec. 7 contains some numerical evaluations of the 
ground state energy. The divergent part is identified as that part of the corre- 
sponding heat kernel expansion which does not vanish for large m (mass of the 
quantum spinor field). After the subtraction of this divergent part the remain- 
ing analytical expression must be transformed in order to lift the regularization. 
A part of the uniform asymptotic expansion of the Jost function is used for this 
procedure. Finally the analytical expression for the ground state energy has been 
evaluated numerically. 

Throughout the paper ^ = c = 1 is used. 

2 The renormalized ground state energy 

— * 

Consider a spinor field on a magnetic background B of the form 

B = ^hir)e- (2.1) 

where $ is constant defining the magnetic fiux, el is a unit vector in the cylindrical 
coordinate system {r,(j),z). The Lagrangian is 

C = -\f^,F^'' + ^ [iY{d, + eA^) - m] ^. (2.2) 

The potential y4(r) 

-> $ air) ^ 

(2.3) 

possesses cylindrical symmetry and the radial part air) is taken to be 

0, r<R^,R2 
<^) = { ^(#-l)^i<^<^2 (2.4) 

1, r>Ri,R2 



3 



The profile function for the magnetic field is 



1 d 



0, r<Ri,R2 

2 



M-) = -^«W= (2-5) 

0, r>R,,R2 

with K = '^^i ^ .The shape of the background can be interpreted geometrically 
as an infinitely long flux tube empty inside. 

We use the zeta-function regularization for the vacuum energy Eq. Since the 
background is static, the following representation holds, 

Eo{s) = -^-n'^J2<"'^ (2.6) 

in) 

where 

^(n) = e(n) = a/^T^ 

are eigenfrequences resp. eigenvalues of energy for one particle states, (the spec- 
trum of the corresponding hamiltonian, see the next section), /z is a parameter 
with the dimension of mass which is introduced to keep the correct dimension of 
energy in this expression. 

For technical reasons we assume that the system is contained in a large but 
finite cylinder of radius R in order to have discrete eigenvalues in the transver- 
sal directions. Because of the translational invariance along the z-axis we can 
separate the 2;-component of momentum, 
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^o(.) = -l^'' J ^ J^ipl + + mr^-. (2.7) 

-oo (") 

Here is it necessary to remark that Eq has the meaning of density per unit length. 
Integrating out Pz we get 

Eo{s) = V^I^4^|^ J2(hnr + m^r-. (2-8) 

In this expression a factor 4 resulting from the summation over spin states and 
over the sign of the one particle energies appeared. The remaining sum over (n) 
is over the eigenvalues k^n) of a two dimensional wave equation for one mode 
in the perpendicular plane. The further transformation of this expression is 



described more detailed in p|, |[T5|. The first step is based on the property 



of logarithmic derivative of eigenfunctions with respect to momentum k. The 
contour integral, 

rik 
Zm 
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is equal to the sum over the spectrum, J^i^n + ^"^Y ^, if the contour 7 encircles 

(n) 

the real positive axis. Then ( |2.8|) becomes 



Ms) = V^-^^l^ E / ^Xk^ + mr-d,lnMkr). (2.9) 

^ 2 ' /=— 00 „ 

7 

Here 0i(r) are the eigenfunctions of hamiltonian H labeled by the orbital mo- 
mentum /. 

In the present formalism instead of wave functions (f)i{r) we use the corre- 
sponding Jost functions fi{k) which contain all necessary informations on the 
spectrum of scattering problem. For the space domain outside of cylindrical 
magnetic background, at r > i?2, the wave functions can be chosen as a linear 
composition of Hankel functions 

Mkr) = ^[Uk)Hl'\kr) + mHl'\kr)], (2.10) 

where fi{k) and fi{k) are the corresponding Jost functions. 

Now after the subtraction of contribution from empty Minkowski space (iden- 
tified as the term divergent ai R —>■ 00, which is independent of the background 
and corresponds to the heat kernel coefficient ao,( p.l5| ) below), and rotation of 
the path towards the imaginary positive axis one obtains |^ 

Eo{s)= -^^1^^ E / dk{k'-m'r^^,lnM^k), (2.11) 

which is a very useful representation of the regularized ground state energy. A 
merit of this representation is the absence of oscillations of the integrand for large 
arguments^. 

To discuss divergences in Eq we need a more general setting. Eq{s) can be 
expressed through the zeta-function of the associated differential operator P 

E,{s) = -^CiP){s~\). (2.12) 
C(-p)(s) admits an integral representation 



00 

1 



Cp{s) = I dt^—-K{t) (2.13) 



^Note that bound states are also taken into account in Eq ( ^.11 ). In the considered problem 
bound states appear if the flux is larger than one flux unit. Strictly speaking, these are zero 
modes located on the lower end of the continuous spectrum, i.e. at k=0 (this is known from 

il)- 
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where the kernel of this integral representation K(t) = **^^(")"^'" ^ is the so 

n 

called heat kernel. The associated "heat kernel expansion" for the kernel K{t) of 



this integral representation at t ^ is ||17 



g -tm 



^ ' n>0 

In accordance to the heat kernel expansion, the coefficients a„ occurred in ( |2.14| ) 



must be for our background [T^ 



aQ = V (2.15) 

is simple an (infinite) volume of configuration space; the part proportional to it 
has been even dropped before during the renormalization procedure as a con- 
tribution of the empty Minkowski space. All other coefficients up to 02 for this 
background must be zero through dimensional reasons and requirement of gauge 
invariance, (for details see the paper 



ai/2 = ai = = 0. (2.16) 
The first nonzero coefficient is known to read 

«2 = ?C (2.17) 
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and only the term with n = 2 contributes to the "divergent" part of energy [p]T[] , 
[0. It will be shown below that only the contribution proportional to 02 in Eq{s) 
contains a simple pole at s ^ and in the case of pure magnetic background is 
proportional to the classical energy 

^dass ^ 1^2 _ 8e^ ^2.18) 
2 3 ^ ' 

In the paper it has been found, that the regularized ground state energy 
of the magnetic fiux tube of finite radius R behaves as — — at i? 00 which 
corresponds to the contribution proportional to 05/2- It holds in accordance to 
the observed fact that in the case of non-smooth backgrounds the heat kernel 
coefficients with half-integer number starting from some value corresponding to 
the smoothness-class of the background are different from zero. It occurs namely 
in our case, where the background potential A{r) is continous and the magnetic 
field h{r) has a discontinuity. The coefficient 05/2 is nonzero in our case, see below 
Sec. 6. The next terms of this asymptotic may be corresponding to 03, and 
a4 delivering terms 

We define the renormalized ground state energy as 

E'-^" = Eo- E^'\ (2.19) 
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where E is obtained from the heat kernel expansion ( p.l4| ) and E^*^"^ fulfills the 



normalization condition 19 



limE'''^" = at m-^oo. (2.20) 



Inserting the heat kernel expansion (p.l4|) into the zeta function and using ( |2.12| ) 
we have 

E^™ = -^(l-2 + /n^), (2.21) 



where 02 is the only non-zero heat kernel coefficient a„, < n < 2, contributing 
to the divergent part of vacuum energy. 

For the renormalized ground state energy the asymptotical dependence on 
powers of m at m — >■ 00 follows from Eq. ( p.l4| ) to have the form 



n>2 

with some coefficients e„. 

In accordance with the interpretation of it must vanish in the limit of 
m — > 00 since it is the energy of vacuum fluctuations. Through the subtraction 
of terms containing all non-negative powers of (which are the terms of heat 
kernel expansion up to 02) the condition (|2.20|) is satisfled automatically. 



Some comments on the subtraction scheme are in order. This scheme is mo- 
tivated by the physical assumption that the quantum fluctuations should vanish 
if the mass of the fluctuating fleld becomes large. The scheme had been used in 



a number of Casimir energy calculations. In |^ it had been shown to be equiv- 
alent to the so called "no tadpole" normalization condition which is common in 
fleld theory. It should be noticed that this scheme does not apply to massless 



fluctuating flelds, for a discussion of this point see ||2T[ . 

So as given by Eq.( |2.19| )is now flnite at s — >■ 0, but it is not possible 
to use this expression immediately since the integral in this limit does not exist; 
i.e. one cannot carry out the analytical continuation to s = 0. In order to 
get a representation where this can be done one can make a trick namely 
add and subtract a part In f^^y^ of the uniform asymptotical expansion for In fi 
in ( p. lip containing a minimal number of terms to provide the convergence of 
the remaining part after the subtraction. Thus we can separate the resulting 
expression into four terms 

E'^"" = {Eo - E^'y"") + {E"' - E*^) (2.23) 

where 



1 ^ 



^asy^ ^2C,y] I dk{k'' - m^y-'^ In fl'^y'^iik) (2.24) 

0rh 
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and redefine two terms in brackets as 



= Eo- E^'y^ (2.25) 

^as ^ ^asym _ ^div ^ ^2.26) 



thus E^ is defined to be 



oo 

^^ = ^ E j dk{e-m') — [lnfi{ik)-\nfrm] (2.27) 



(=— oo 



The "finite" part E^^ of the renormahzed vacuum energy is not only well defined 
at s ^ 0, but also provides a representation well suited for numerical analysis. 
The analytical continuation for i?"* at this limit will be constructed below. 

3 Solution of the Dirac equation 

We consider a spinor quantum field ip in the background of the classical magnetic 
flux. We start with the Dirac equation for this fleld 



^r^-m-eYA,\^P = (3.1) 



with the electromagnetic potential ( [^.3| ). The ip in ( ^.1| ) is a spinor quantum field 
of mass m and charge e, interacting with A ( p.3|) . The gamma matrices in our 
representation are chosen to be the same as in 



o"3 1 1 _ / 1 2 _ / —io-i 1 3 _ / 1 



^ ^ -aa J ^ \ -ta2 j ^ \ la^ j ^ {-10 

(3.2) 

Now we follow the standard procedure and separate the variables. Using the 
ansatz 

0(r, if, z) = e-'P«"°e-*^'^"^(r, ip) (3.3) 



we obtain the equation for the 4- component spinor \E' 

P3(^3 Po + L + mas 



X 



^ = 0, (3.4) 



2 

where L = i ai{di + ieAi), 0, x are the two-component spinors. 

i=l 
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Respecting the translational invariance of the system it is sufficient to solve 
the equations only for ^3 = 0. Then one of the two decoupled equations (0, x) 
reads 

' po — m 



d_ _ l-Sajr) 
dr r 







d_ _ l+l-Sa{r) 
dr r 



(3.5) 



Po + m 



(here the standard ansatz 

V''(r)e-^''^ 

has been used, / is the orbital quantum number), the radial part denoted as 



(3.6) 



^(r) 




and 



2^' 



(3.7) 



(3.^ 



In order to use later the symmetry properties of the Jost function we redefine 
the parameter / in (|3.6| ) (orbital number) as u according to 



l + \ f(yrl = 0,1,2,... 



V 



-/ - i for I 



(3.9) 



-2,... 



For the given construction of potential ( p.4|) we get three equations and three 
types of solutions for 3 areas of space respectively 

Domain I: r < Ri. The free wave equation (the magnetic flux is zero), the 
solutions are the Bessel functions Jj^_^i{kr) 



(r) 



^r(r) 



<|.+(r) 



1 



y/po - mJ^^i{kr) 



Po = y/m"^ — k"^ 

This solution of equation ( p.5| ) in the domain I is chosen to be the so called regular 
solution which is defined as to coincide for r — * with the free solution. 
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Domain II: Ri < r < R2 The equation with a homogeneous magnetic field 
has the solution 



Il.r 



4 



Il.r 



4 



c: 



(po+>»)».a+ig-^ F Q _ hi a + 2-^) 

2(5+1)' l-^Ul 2/3'"^^' 2 ) 



r"e-^iF,(-|,ci + l;^) 



;^)-^e-V,Fi(-§-a,l-a;^^ 



(3.10) 



$/j(r) is the same, but a replaced by a 
here: q; = z/ — i + 4for/>0 



- z/ - i for / < 



iFi is a confiuent hypergeometric function |p5 |. 

The coefficients Cj, Cr are some constants that will be irrelevant for expressing 
the Jost function. The indices u and 1 denote "upper" and "lower" components 
of spinor respectively. The lower index "i" or "r" corresponds to "regular" or 
"irregular" part of solutions dependent on the behaviour of the function if con- 
tinued to r = 0. If the external background vanishes {6 0), contributions of 
irregular parts to the solution disappear. 

Domain III: r > R2 The free wave equation outside of the magnetic fiux has the 
solutions 



^iiiir) 



2 y/po - m 



rniAr) 



iiiA 



2 y/po - m 



rniAr) 



^JAik) 



Kl^^sikr) 



+ ^L-ik) 



H^^\ Jkr) 
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2 VPo - vfi 



r) 



1+ 



L •^m.r-v 



+ \ft{k) ^ 



and fu{k) is the desired Jost function in accordance to the definition ( |2.10|) . The 
functions fi{—k) and fi{k) are conjugated to each other because of the choice of 
the solution in the domain I to be regular (Rl). 



4 The exact and asymptotic Jost function for 
the background 

To obtain a Jost function we need to impose certain matching conditions for the 
spinor wave at the boundaries between I-II and II-III domains. It is known that 
for a continuous potential a(r) as we consider here, the spinor wave function must 
be continuous, hence for its components it holds 



r=Ri 



r=R2 



Resolving these equations for we obtain, that 



(4.1) 



A 
B 



o A + B 
WniWn 



Rl 



Rl 



l± J.U± 



Il.irill.r) 



R2 



R2 



(4.2) 



The denominator of this expression can be written using the Wronskians of 
hypergeometric and Bessel functions as follows 



Win = Wn..4ii.r - 4ii.^m.r) 



,Po + m.i Ai 

2 



R2 



Pq — m 7rA;i?2 



-OL-l 



Rl 



Po 



m 



(4.3) 
(4.4) 
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To calculate Eq we need this function with imaginary argument ( p.ll| ) so we 
need to replace k by ik. As a result we obtain the new expression for f{ik) that 
contains now modified Bessel functions I^^i 

functions K,,,i^: instead of H^^''^i 



instead of J,,_i_i and modified Bessel 

2 '•^^2 



nkRiR2e-^/''- 



iirS 
2 



X 



kRi 



2i> + 1 ^ 2/3' 2 



3 pRl 



2i> + l 1 A;2 1 

^ ^^^^2"''+2^' 2 



K.sikR2)- 



k iFii z/ H , u; — 

^ ^^2 2/3' 2 '2 



(A;i?2) 



kR 



P . 3 pRl 



2^ ^^^^' + 2^' ^"+2 



^^^^2^' '+2 



1 /3i?2 



where u = u + -^73^, and 
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-2 + ^ 

29-1 ^ 2 2f3 2 2 ' ^^+5^ ^ 

1 P 3 



(4.5) 



X 



'2u-l e 3 _ pRl 

k iFi{—-, z/; 



is:.. 



k^ 



1 



K,^i.sikR2) + 



I k'^ 3 /3i?i, „ , 



X 



where i' 



- and 

Now we need to obtain the asymptotic Jost function and asymptotic part of 
energy. We use the representation of the asymptotic Jost function In fas{ik) in 
the following form (see App.) 



n=i 7=1 



(4.6) 



where t = (1 + and Xnj ( given explicitly in App. A) are represented in 

terms of r, 6, a(r) and their derivatives. 

This expansion had been obtained in [T^ by iterations of Lippmann-Schwinger 
equation up to the order u^^. In general it is possible to obtain higher orders 
using this formalism. However as the calculations |TB[ showed, the complication of 
the involved expressions increases very fast. It is remarkable that this expression 
does not contain a term with power z/"^. In the finite part of the energy the 
corresponding term is canceled in the sum of terms corresponding to positive and 
negative orbital momenta z/ as well. The absence of the power z/"^ is a succession 
of the zero heat kernel coefficient 03/2. Also it is a non-trivial fact that both, the 
fourth power of the magnetic flux 6 and the second one is present. 
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It can be checked numerically that Inf""^ ( |4.6| ) is indeed the uniform asymp- 
totic expansion of logarithm of ( [4.6|) for k, u- large, k/u = 2;-fixed. 

5 The finite and asymptotic parts of vacuum 
energy 

The expression ( [4.6|) for the uniform asymptotic of Jost function, substituted into 
the expression for E""^ ( p.26| ) yields 

OO OO 2 g OO 

E'- = 2C, Y: jik(k^.,n-Y-'^J'tY^Y.j^"'^,-E'- (5.1) 







i>=l/2,3/2,...;; n=l i=l 

with the constant = ^[1 + s(-l + 2Zri2/i)]. E''™ as defined in (g]^ is 



^ = a rp2 p2, (- - 2 + In ^) (5.2) 
67r(_K2 — iti) s 

with the coefficient 02 according to (|2.17|) has been used in the form 



00 

2^2 8^x2 /"^ ^2 167r52 



a2 = = —5^ / drrhirV = . ^ (5.3) 



Then the sum over u can be now transformed into two integrals by means 
of Abel-Plana formula ( |9.1|) . The first one cancels the E'^™ exactly. Then we 
integrate the second one over k using the identities (|9.2| , |97^) . It gives the form 



n=l j=n Q 



rm) (5.4) 



with 



T{s + 3/2-1) -I 
dry ( iiuy-'^ (-iz/)^"" 



1 + eXp{2T^v) V [1 + (iz//x)2]^+i/2-l [1 + (_iiy/a;)2]«+i/2-l 

u 

In order to obtain an analytical continuation in s = of each term of the sum 
we integrate it over v by parts several times till the divergency at z/ = through 
the power s + j/2 — 1 abrogates, and after that we can perform the integration 
over r. Further we integrate by parts resulting in the relations raa' —\a^rdr 
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and r'^aa" — a'^ + ^a^rd^r (that hold because of the continuity of the potential 
a(r)), and we obtain finally the form 



—4 r fir 

= _I / ^[S^a{rygi{rm) - 6\^a'{rfg2{rm) + (5'a(r)^^3(rm)] (5.6) 





oo 



with gi{x) = J dv\/v'^ — x^jiiy) that will be calculated numerically (see the 



next section), the functions fiiv) are shown explicitly in 

The finite part of the ground state energy can be finally represented as 
follows 

oo 

^^ = ^ E f dk{e~m')^[lnf;^itk) + lnf-{tk)-2\nritk)] (5.7) 

^ u=l/2,3/2,...i 

In order to use this form for numerical evaluations we integrate by parts and get 



oo 



{5.t 



i/=l/2,3/2,...; 

(the logarithmic expression in square brackets in (|5.7|) denoted as Inf^^^) 



6 Higher orders of the uniform asymptotic ex- 
pansion of the Jost function and the heat ker- 
nel coefficient 

The background considered in this paper has singular surfaces where the magnetic 
field jumps. The heat kernel expansion for the case of singularities concentrated 
at surfaces has been considered in [|^|, |TB[, [^]. Although the general 
analysis of |18| is valid for our background, an explicit expression for 05/2 has not 
been calculated yet. 

We can use our obtained Jost function ( [4.6| ) to calculate the coefficient 05/2 
in the heat kernel expansion ( p. 141) . 

Suppose we have obtained the value of Eq{s) ( p.llD in the point s = —1/2. It 
follows from (gl2| - |]l|) that 

Eo{s 



m 



2(47r)3/2 r(s 



2n 



T(s -2 + n) 



(6.1) 



and at the limit of s 
in the sum. 



2^ n=0,i,l,... 

— I only the term containing 05/2 remains to be nonzero 



/i 



2s 



2(47r)3/2 
15 



«5/2 



(6.2) 



From the other hand we have ( |2.11| , |2.24| ) 

oo 

Eo{s) = Cs Yl I - ^')'''^/- W = ^Ms)- (6.3) 

i/=l/2,3/2,...;; 

Here we substitute the exact Jost function by its uniform asymptotic represented 
in the form 



n=l,3,4,5... 

where the coefficients hn are functions of ti = (1 + {^^Y)^^^"^ , ^2 = (1 + 
^fci?2 ^2^-1/2^ the power z/~^ is absent, as noticed above 



_ /i^T(s-l) -4sin(7r3) 



At the hmit s ^ — it yields 



2' Stt 

and therefore we obtain for 05/2 



1 4 

Eo{--) = —Res^^_ih{s) (6.6) 



^5/2 



''^^i^es.^iMs). (6.7) 



3 

To obtain the exphcit form of ReSg^_ih{s) we use the uniform asymptotic 



expansion of the Jost function ( p.8|) . The terms /in(^i;^2) in |6.4] can be obtained 



either by iterations of Lippmann-Schwinger equation (see ||T3[, §]) or by using 



the exphcit form of the Jost function as well. All the further terms up from 
/i4 are produced from the explicit form of the Jost function (^4.6|) because of 
complication of the first way for higher orders n (see the remark to in Sec. 6). 
Namely, we obtain several higher orders 1/z/ of uniform asymptotic expansion 
for special functions 1^, Ky^i Fi which the exact Jost function ( [4.6| ) consists of (it 
can be done starting with the explicit form for two first orders and executing the 
recursive algorithm several times |2^), then after substitution of each of functions 
/, K,i Fi by its corresponding uniform asymptotic expansion and separation of 
powers of u we arrive at the form ( |9.8|) . The coefficients hn{ti,t2),n = 1, 4 are 
given in the Appendix (|978|). 

If the function /i„(ti,t2) is a polynomial over ti,t2, so we can consider some 
term P,(t = ti,i = 1,2) of it. Notice, that for /ii(ti, ^2), ^3(^1, ^2) it is not the 
case, but we can treat the terms of kind and ' ^-integer, as an infinite 
sum of powers t. 

Apart from the construction of t = ti,t2 (3:4) they are strictly positive and 

00 

less than 1, therefore the series = Yl^^ converges regular and uniform and 

1=0 
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respecting that we have only finite integer powers k so that [ytj] = ( X] ^* ) 

converges as well, the sum over i can be interchanged with the one over v (Prince- 
heim's Theorem) thus the following procedure is valid. Performing the sum 



h{s) 



oo 

1/2.3/2.... 1, 71=1,3,4,5 



(6. 



i/=l/2,3/2,...; 

by meaning of (|9.1| ) we obtain the sum of two parts 
h{s) 



oo oo 



(6.9) 



dv 



where i = ti^2- The first summand of ( |6.10| ) gives for each power j of t = {^1,^2} 
(using the (fy)) 



m' 



r(2-s)r(i±p)r(s + 



n-3 ' 
2 ' 



2 {RmY-^V{]/2) ^^'^^^ 

(where R denotes Ri and R2 respectively) at the limit s — > — | only the terms 
corresponding to n = 2 and n = 4 (which will be calculated below) have a pole. 
For the second part of the Abel-Plana formula ( p.lO| ) we have 



dv r(2-s)r(s + j/2-l) 



1 + e 



2'ITU 



ji?4-2.r(_^72) 



X [(zz/)2+^-"[(mi?)2 + {ivfY-^/^-' 



-IV 



,2+j-n 



[{niRf 



-IV 



(6.11) 

|2jl-i/2-sj 



It can be seen using the Taylor expansions at s — >■ — | of Gamma functions 
entering the ( |6.12 ), that only the terms containing -1, 1, and 3 powers of ti can 
contribute to the residuum at s ^ — But the expression in square brackets in 
( |6.12| ) can be performed as 



-z/2(z/2 - {mrf)[{ivy-''e'^^^-^/^-'^ - (-iz/)^-"e-*"(^~^'/2~")] 



(6.12) 



thus for j —n even the expression in square brackets produces (dropping the non 
sufficient coefficient ±2 or ±2z): 

for j = — 1, 1 : sin7r(l — j/2 — s) = cosTTS 

for j = 3 : sin7r(l — j/2 — s) = — cosTTS 



and for j — n odd respectively 
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for j = ±1 : cos 7r(l — j/2 — s) = ± sin TTS 
for j = 3 : cos7r(l — j/2 — s) = — sinvrs 



But for s ^ — I these functions behave as 



^'1 I „^2 



siUTTS ~ — 1 + + 

COSTTS ~ 7r(| + s) 

and it means that only the contribution of terms with odd powers j — noi iv could 
survive and these corresponds for the possible values of j to z/~^, .... In 

fact the coefficient /i2(s) at is zero, and all other possible terms up from z/~^ 
does not contain any powers of ti lower than 4; one can see it for example in the 
explicit form of the uniform asymptotical expansion of special functions entering 
the ln/t,(zA;) (|4.6|). Therefore the second summand of ( |6.10| ) does not produce 
any contribution to the ReSg_,__ih{s) Thus we have that only the contribution 
from the first summand of (|6.1CI|) remains, and since the term of = 2 is zero the 
searched residuum resulting from the term of n = 4 is: 

Res.__J.is) = + fi.) E (6.13) 



where the ( |6.10| ) and the explicit form of /i4(ti, ^2) = j(j^T^:]^[Rii'tt~'ti) + R^i'^t' 



Res,^_.h{s) = — (6-14) 



^2)] have been used. Finally we have 

157r 6^{Ri + R2 
128 {Rl - RlY 

and therefore 

h7r'l H\R^ + R,) 
"^/^ = {Rl - RIY ^^-'^^ 

This is the heat kernel coefficient 05/2 for the configuration of the magnetic back- 
ground field as given by Eqs. (p73| - |2.5| ). 



We can calculate the heat kernel coefficient 05/2 in a more general situation 
when the magnetic field jumps on an arbitrary surface S. The coeffcients a„ for 
n = 1/2, 2 can be read off rather general expressions of the paper [|l^. Let 



be values of the magnetic field on two sides of S. According the analysis of [0, 
the coeffcient must be an integral over S of a local invarinat of canonical 
mass dimension 4, which is symmetric under the exchange of and B~ and 
which vanishes if B~^ = B~ (i.e. when the singularity disappears). There is only 
one such invariant which gives rise to the following expresssion: 

05/2 = e j {B+ - B-fdii{^), (6.16) 

where the integration goes over the surface S and B"^ are the values of the 
magnetic field on both sides of S in the given point. The yet undefined constant 
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^ can be found using Eq.( |6.15| ) which constitutes a special case of ( [6.161 ). Here 
the surface S consists of two circles in the (X, y)-plane so that 



(fi+ - B-yd^iX) = 2tt{Ri + R2)B' 



(6.17) 



where the jump (5+ — B ) is just the value of i? at r G [Ri,R2]. With this, 
Eq. (|6.16|) takes the form 



a5/2 = 27re2e^'(^i + ^2). 
On the other side, from Eg. ( p.l|) we have 



$ = 5 j Bd^x = nB{Rl - Rl) 

relRuR2] 

and with Eg. ( |3.8D from Eg. ( |6.15[ ) it follows that 

5 



0-5/2 



__^2^2 R1+R2 



8 A im-Ri?' 

comparing (|6.18|) with ( |6.20|) we get 

5 



167r3/2- 



(6.18) 



(6.19) 



(6.20) 



(6.21) 



7 Graphics and numerical results 

The asymptotic part of the ground state energy as given by Egs. ( [5. 4] ) and ( |5.6| ) 
can be represented as sum of a part proportional to the second power of the flux 
and one proportional to the fourth power, 

= --[ei{Ri,R2)6^ + e2{Ri,R2)S% (7.1) 

TT 

The corresponding coefficients ei(i?i,i?2) and e2{Ri, R2) can be calculated nu- 
merically without problems. They are shown as functions of R2 for fixed Ri in 
Fig.! (multiplied by Rj). 

The finite part of the ground state energy, Ef, is used in the form as given 
by Eg. ( p.8|) . Let Ejiy^k) denote the function to be integrated and summed 
over in that expression. In Fig.^ it is shown as function of k for several values 
of the orbital momentum v. In order the make the behavior better visible it is 
multiplied there by v'^k^. These functions are smooth for all values of fc, starting 
from some finite values at A; = 0. For large /c, the functions v^k'^Ef{v^ k) shown 
here tend to a constant thus the integral over k is convergent. All integrals 
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Figure 1: The asymptotical part of energy as function of the outer radius for 
= 0.00001 resp. Ri = 1 

have been truncated at k = 1500. The error caused by this is quite small and 
does not change the results shown in the Table 1. These integrals we denote by 
Ef{u). They are shown as function of v in Fig.^ in a logarithmic scale. Again, we 
multiplied by a power of argument, here by v"^, in order to make the behavior for 
large v visible. It is seen that v'^Ef{v) tends to a constant so that the sum over v 
is convergent. The sum is taken up to z/ = 232.5 and again the remainder is small. 
The calculations have been performed for several values of the parameters. The 




-1 . 5 



Figure 2: The contribution Ef{v,k) multiplied by v^k^ of the individual radial 
momenta for several orbital momenta to the finite part of the ground state energy 
for Ri = 0.0001, R2 = 1,6 = 3. 

results are displayed in Table. 1. The computations are performed with an adapted 
arithmetical precision. In intermediate steps compensations between sometimes 
very large quantities appeared. The precision was adapted accordingly. For 
example, for Ri = 0.99, u = 250.5, k = 1600 as much as 1404 decimal positions 
have been necessary to get at least 16 digits precision of the integrand Ef{v, k). 
This was a factor causing large computation time. 
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v'E^ (v) 

0.6- 
0.4 ■ 
0.2 ■ 




Figure 3: The contribution Ef{y) of the individual orbital momenta to the finite 
part of the ground state energy multiplied by for — 0.0001,i?2 = 1 and 
(5 = 3. 

Table 1. The numerical evaluations for several values of i?i and 5 



i?i=0.0001 



5 








j^class 




0.5 


-0.0130152 


0.00152886 


-0.0114863 


3.14159 


3.13011 


1. 


-0.052403 


0.00509591 


-0.0473071 


12.5664 


12.5191 


3. 


-0.450093 


-0.052012 


-0.502105 


113.097 


112.595 


6. 


-0.91953 


-1.52936 


-2.44889 


452.389 


449.94 


10. 


4.85954 


-12.9483 


-8.08871 


1256.64 


1248.55 


15. 


46.3333 


-67.3661 


-21.0328 


2827.43 


2806.4 


21. 


215.273 


-261.526 


-46.2534 


5541.77 


5495.52 


30. 


987.62 


-1095.29 


-107.666 


11309.7 


11202.1 



Ri=0.9 



6 






'^ren 


j^class 


^tot 


0.5 


-0.278679 


-0.00105651 


-0.279736 


16.5347 


16.255 


1. 


-1.11431 


-0.00451479 


-1.11882 


66.1388 


65.02 


3. 


-10.0633 


-0.0683525 


-10.1317 


595.249 


585.117 


6. 


-40.4372 


-0.647622 


-41.0848 


2381. 


2339.92 


10. 


-112.17 


-4.2629 


-116.433 


6613.88 


6497.45 


20. 


-423.279 


-63.2506 


-486.53 


26455.5 


25969. 


40. 


-1079.05 


-992.187 


-2071.21 


105822. 


103751. 
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i?i=0.95 



6 








j^class 


J^tot 


0.5 


-0.72166 


-0.00212769 


-0.723788 


32.2215 


31.4977 


1. 


-2.86633 


-0.00877711 


-2.87511 


128.886 


126.011 


3. 


-25.8293 


-0.104564 


-25.9339 


1159.97 


1134.04 


10. 


-287.834 


-4.39364 


-292.228 


12888.6 


12596.4 



i?i=0.99 



5 








j^class 


J^tot 


0.5 


-5.59211 


-0.0106385 


-5.60275 


157.869 


152.266 


1. 


-23.2975 


-0.0428038 


-23.3403 


631.476 


608.136 


3. 


-202.057 


-0.409235 


-202.466 


5683.28 


5480.82 



/?i=0.997 



5 


E^ 






j^class 


J^tot 


1. 


-91.996 


-0.142036 


-92.138 


2097.54 


2005.4 



i?i =0.999 



s 


Ef 






j^class 


J^tot 


1. 


-311.182 


-0.425555 


-311.608 


6286.33 


5974.72 



8 Conclusions and Discussions 

In the preceding sections the ground state energy for a spinor in the background 
of a rectangular shaped flux tube had been numerically calculated. The cor- 
responding Jost function had been written down explicitly, also its asymptotic 
part. The numerical calculation required work with high arithmetic precision. 
The results are displayed mainly in Table 1. For small inner radius of the flux 
the results are close to them of where the same problem for a flux tube with 
homogeneous magnetic field inside, which corresponds to i?i = here, was con- 
sidered. Especially, it is seen that for large flux 5 there is a compensation of the 
^'^-contribution between the finite and the asymptotic parts of the ground state 
energy leaving a behaviour proportional to 5"^ In 5 as shown in Fig.§. Here the 
asymptotic part gave an essential contribution. The ground state energy remains 
negative, but numerically small. Only for very large flux it could overturn the 
corresponding classical energy, but these values of the flux are clearly unphysical. 

For values of the inner radius close to the outer one, Ri ^ 1 (where we have 
put i?2 = 1) the picture changes. Here the vacuum energy grows faster than the 
classical one. Generally, both diverge proportional to (1 — -Ri)^^, the classical 
energy is equal to ^^ciass = 5^27r(l — -Ri)~^. The vacuum energy, multiplied by 
(1 — is shown in Fig.|^ in a logarithmic scale. It is negative and growing a 
bit faster than the classical one which would be a constant in this plot. Here the 
asymptotic part of the ground state energy becomes increasingly unimportant 
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Figure 4: The ground state energy divided by S"^ In S as function of S. 

(see Table 1). The question whether the vacuum energy for sufficiently small 
(1 — Ri) may become larger than the classical one cannot be answered by the 
numerical results obtained. The problem is that the computations become very 
time consuming because of the increasing precision which is required. Also, one 
has to take higher k and u into account. The weakening of the growth for Ri = 
0.997 and Ri = 0.999 seen in Fig.^] may be caused just by this reason. Here one 
has to note that the integrand is for large k and u always negative (see Fig.^) so 
that dropping some part (as we did within the numerical procedure) diminishes 
the result. So, as a result, we cannot exclude from the given calculation that 
the vacuum energy grows for a strong background field faster than the classical 
energy. 

Further work is necessary to better understand these questions. An improve- 
ment of the numerical procedure is certainly desirable. It could go along two 
lines. First, in the calculation of the Jost function the compensation of large ex- 
ponentials should be avoided by taking them into account analytically. Second, in 
the compensation between the logarithm of the Jost function and its asymptotic 
expansion in the integrand of Ef in Eq. (|5.7|j5^ ) higher orders of the asymptotic 
expansion could be used. However, for this reason one would have to continue 
the procedure invented in |jl3[ for this expansion using the Lippman-Schwinger 



equation or to find some equivalent procedure. 
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Figure 5: The vacuum energy multiplied by (1 — -Ri) in a logarithmic scale for 
5 = 1. 
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9 Appendix 

The sum over v has been transformed to integrals using the Abel-Plana formula 
as follows: 

oo oo 

f(, W -^'-> (9.1) 

'=0 

The integrations over u and k can be done using identities: 



oo oo 



m 



* ' « " r(i) (1 + 

The expansion in powers of v for logarithm of asymptotic Jost function can be 
obtained in the form (see0) 
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f dr 

n=l j=l 



dr „ P 



(9.4) 



where 

Xi,i = 


-Xi,3 = X2,6 = |(a5)2 


-^2,2 = 




X2,4 = 


i52(_3a2 + raa') 


-^3,3 = 


i52(a2 _ ^^^/ + 1^2^^// _ i^2^4) 


-^3,5 = 


|52(_|^2 ^ 7^^^/ _ ^2^^// ^ g^2^4) 


-^3,7 = 


|52(35a2-5raa'-552^^) 


-^3,9 = 


Id 



For the representation of Eas as 5.6 



oo 

-4 /"(ir, 



E"' = j '^[5'^a{rfgi{rm) - 6\^a{rfg2{rm) + 5^a(r)^^3(rm)] (9.5) 



here fi are 

fiix) = - ^/i,3(a;) + ^/3,3 - ^/3,5(a;) + y/3,7(a;) - ^/3,9(a;) 

11 7 5 

- 2^5x(--/3,3(a;) + -/3,5(a;) - -/3,7(a;)) 

+ ^a;9^(|/3,3(a;) - ^/3,5(x)), 

/2(X) = i(/3,3(x)-/3,5(x)), 

fsix) = -^(/3,3(x)-6/3,5(x)+5/3,7(x)). (9.6) 

with fi,j are 
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/l,3(x) 



1 + e^^^^ 

1 1 w 



1 + e 



11 X 

hri^) = 4(-(-(t^^)')')' 

15 X X 1 + e^'^^ 

105 X X X 1 + e^^^ 

The asymtotic of the logarithmic Jost function can be obtained in the form: 
In fi,{ik) = Y,hn{ti,t2),{see |6]|) with: 

n 

hi = ^\^[RiHi{2 + ti)(l + - 3Ri^R2\l + ti)%{l + ta) + i?2^(l + ti)\{l + 2t2)] 



h, = Ax^[2R,Hi\l+t2f + Ri'R2\l + ti)%'{l + t2)\3t2'-4)- 

D 

R2\l + ti)W3t2^ + 6t2^ - t2^ - 8t2 - 2)] + 

^\^[ARihi%A + ti)(l + ta)' - 5Ri'R2^{l + ti)%%l + ta)' + 

15i?/i?2^(l + tl)'t2'(l + t2)'(t2' + 2t2 - 1) - 

5Ri^R2%l + ti)%^{3t2^ + I2t2^ + 6t2^ - At2 - 1) + 

i?2^(l + tl)'t2^(5t2^ + 20t2^ - 4t2 - 1)] 

, _ 5^[RiHi\l-ti^) + R2%\l-t2^)] 

iIa — T\ 

mi - Rif 

^ = (Rl-Rl)(l+ti)il+t2) ^^-^^ 
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